% retorna una con ceros debajo de la diagonal
% en este caso con los multiplicadores debajo de la diagonal.
function res_matriz = triangular_matriz(matriz)
    % calcula el limite de la iteracion, por si la matriz no es cuadrada
    [fil,col] = size (matriz);  
    limite = min (fil,col);
    for j = 1:limite %iteracion por columnas
        display(datestr(clock));
        for i = j+1:fil %iteracion por filas
            % si ya hay un cero en la posicion evaluada
            % debe pasar a la siguiente iteracion
            % calcula el multiplicador 'm'
            m = matriz (i,j)/matriz (j,j);
            % guarda el multiplicador
            %matriz (i,j) = m; 
            % termina de operar con la fila
            for k = j+1:col
                if  abs(matriz(j,k))> 0
                    matriz (i,k) = matriz (i,k) - (m * matriz(j,k));
                end
            end            
        end
    end
    res_matriz = matriz;
end
    